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We provide a systematic analysis of the physical generation of single- and two-qubit quantum 
operations from Hamiltonians available in various quantum systems for scalable quantum infor- 
mation processing. We show that generation of one-qubit operations can be transformed into a 
steering problem on the Bloch sphere, whereas the two-qubit problem can be generally transformed 
into a steering problem in a tetrahedron representing all the local equivalence classes of two-qubit 
operations (the Weyl chamber). We use this approach to investigate several physical examples for 

- - - the generation of two-qubit operations. The steering approach provides useful guidance for the 

'nI realization of various quantum computation schemes. 
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INTRODUCTION 



Physical implementation of quantum information processing and quantum computation usually begins with a cou- 
pled quantum mechanical system and requires that the available Hamiltonian for that system be controlled to generate 
C^l , desired quantum operations from which the quantum algorithm is constructed ^Ij, |2|j . It is therefore a fundamental 

(N 

issue to produce quantum operations from the Hamiltonian that is provided by or accessible to the physical quantum 
^ , system. In this work we formulate the generation of these operations as a problem in control theory and show that 
^ ■ this language is beneficial to developing efficient physical implementation of quantum operations for computation and 
information processing. 

\^ The postulates of quantum mechanics assert that the state of a quantum system is completely described, at time 

T-H ' i, by a unit vector \ip{t)) in a Hilbert space 3]. The evolution of |V'(i)) is determined by = ^(OIV'(O)): 

, where U{t) is a unitary evolution operator (or propagator). The dynamics of U{t) is given by Schrodinger's equation 
ihU{t) = H{v)U{t) with J7(0) = /, where H is the Hamiltonian of the system, and v the external control field. 
The generation of quantum operation means to find a control v such that the trajectory generated by Schrodinger's 
■ equation can achieve a prescribed target unitary operator Ut at certain final time. This can be posed as a steering 
\ problem in control theory. The state space of this control system is the unitary Lie group U(2"'), which consists of 
all the quantum operations on n-qubit system, and the dynamics of quantum operation described by Schrodinger's 
\ equation is a right invariant vector field on the unitary Lie group U(2"). 
^ • Steering on Lie groups has been studied for a long time, and it still remains a very difficult problem in control 
^ ■ theory Q. It is well-known that when the control distribution is bracket generating, the system is controllable, i.e., 
Q^, there exists a control law that steers the system from any initial state to any final state in the state space ^ 0| . 
^ Nevertheless, this result does not provide a constructive way to find the control law. One useful mathematical tool 
• T—i to solve steering problem is to use Pontryagin's maximum principle, which gives a set of differential equations that 
' have to be satisfied by the optimal trajectory [tIIM lol Il0| . Quite often however, this approach leads to a two-point 
^ \ fixed-boundary problem which has no analytic solutions . 
■ - - 1 In this paper, we develop a general approach to design the control law for generation of an arbitrary quantum 
operation from the available Hamiltonians in quantum mechanical systems. Because any quantum operation on any 
arbitrary large quantum system can be decomposed as a combination of single- and two-qubit quantum operations jj] , 
we only need to study the control of single two-level systems and pairs of coupled two-level systems, which amounts 
to the steering problem on the unitary Lie groups U(2) and U(4) respectively. The key idea of the approach here 
is based on the observation that in many real quantum systems there exists a certain subset of quantum operations 
that can be readily generated. This observation leads naturally to the notion of control on a quotient space, that is, 
instead of controlling the original high dimensional complicated quantum system, we study a reduced control problem 
on a quotient space obtained from an equivalence relation that is defined for two quantum operations if they differ 
only by quantum operations in that specific subset. Therefore, to generate a desired target quantum operation, we 
can equivalently generate a quantum operation that differs from the target by a quantum operation in the easily 
achieved subset. One key issue in this approach is to find the structure of the quotient space and the reduced control 
trajectory on this quotient space. This makes it possible to study the generation of quantum operations by solving 
a steering problem on the quotient space. We present solutions here for steering on the quotient space of both U(2) 
and of U(4). The solutions for U(2) are related to those known from the NMR literature, but are developed here in 
a unified steering framework that also encompasses the solutions for the more complex U(2). 

We start with the generation of single-qubit operations, which can be viewed as a control problem on a two-level 
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system. We know that all the single-qubit operations form the unitary Lie group U(2), and that any single-qubit 
operation can be described as a rotation on the Bloch sphere of a certain angle about a specified axis. In many 
two-level quantum systems, it is often easy to generate a rotation operator Rn about a fixed axis n by appropriate 
choice of a static control field. By making a change of coordinates one can transform Rn into i?^, i.e., the rotation 
operator about the z-axis. Consequently, to implement an arbitrary single-qubit quantum operation Ut, we only 
need to implement a quantum operation Ui that differs from Ut by a z rotation operator from the right, that is, 
Ui = UtRz{0), with 9 g [0, 27r). This defines an equivalence relation on all the single-qubit quantum operations, and 
the quotient space obtained from this equivalence relation is nothing but the Bloch sphere. Therefore, the generation 
of single-qubit operations can be reduced to a steering problem on the Bloch sphere. For the case when the control is 
an oscillating field perpendicular to the static field, we derive the control laws for generation of arbitrary single-qubit 
operations that are familiar from the NMR literature. We also prove that the system is controllable if and only 
if the oscillating field is tuned at the resonant frequency. We then study the case when the oscillating field is not 
perpendicular to the static field. For both of these two cases, we show that when at resonant frequency, the total 
time to implement a given operation is quantized. This is a particularly important issue for the implementation of 
local unitary operations in coupled two-qubit systems, for which two single-qubit operations must be generated at the 
same time. We also show that when the control variable is varied within a certain range, an efficient way to generate 
any single-qubit operation is to simply alternate between the two extremal values of the control parameter. A control 
law design based on the Bloch sphere steering technique thus provides a complete solution to the generation of any 
arbitrary local unitary operations from a given Hamiltonian. 

We then consider the generation of two-qubit quantum operations. Physically, this is a control problem on two 
coupled two-level systems. It is much more difficult to control a coupled qubit system to generate the desired two-qubit 
quantum operation, because of the higher dimensionality and more complicated Lie group structure. In Ref. |1C| . 
Khaneja et al. derived the time optimal solution when the Hamiltonian is purely nonlocal. However, when the 
Hamiltonian contains both local and nonlocal terms with c omp arable magnitudes, there are no systematic methods 
to derive the control laws other than numerical procedures [13, ^tii^ paper, we provide a systematic approach 

to solve this problem by using the notion of control on a quotient space. Our control strategy employs two steps to 
generate an arbitrary two-qubit operation Ut- First, we find an efficient control law to generate a nonlocal quantum 
operation Ui that differs from Ut only by local unitary operations fci and ^2, where each of these is the tensor 
product of two single-qubit quantum operations. Then we consider generation of the local operations fci and k2 that 
is required in order to arrive at Ut in the computational basis. This will often require that the single-qubit operations 
be implemented simultaneously (see above). 

A basic principle in this control design strategy is that we implement the local and two-qubit operations sepa- 
rately. This is often feasible in many physical applications, where the single-qubit operations are generated by local 
Hamiltonians that can be independently controlled, e.g., by external fields, so that we can readily implement any 
local unitaryoperation. This separation leads naturally to the concept of local equivalence on all two-qubit quantum 
operations In Ref. we derived the geometric representation of this equivalence relation as a tetrahedron, 

referred to as the Weyl chamber. Each point in the tetrahedron represents a local equivalence class of some nonlocal 
two-qubit operations. The generation of a two-qubit quantum operation can then be treated as a steering problem 
in this tetrahedron. We have shown some examples of this steering technique 15], when the Hamiltonian is given 
in some special cases, e.g., purely nonlocal. In this paper, we will exploit this idea to more systematically study the 
generation of two-qubit operations for a variety of Hamiltonians that include both nonlocal and local terms and that 
derive from physical systems relevant to implementation of quantum information processing, in particular to solid 
state implementations. 



II. GENERATION OF SINGLE-QUBIT OPERATIONS: STEERING ON BLOCH SPHERE 

In this section, we study the generation of single-qubit operations as a control problem of unitary transformations 
on two-level quantum systems. The time evolution operator for a general two-level system can be written as follows: 

lU ^{Hd + vH,)U, UiO)^I, (1) 

where Hd is the drift Hamiltonian, He the control Hamiltonian, and v the control. Here and He are both Hermitian 
matrices on C^^^. Upon neglecting the global phase term, all single-qubit quantum operations form a Lie group SU(2): 

SU(2) = {Ue C^^2 . jjui ^ J ^i^^ u ^ ^ 

Our task is to find the control v that will drive the system from the initial operation U{Q) = / to a prescribed 
target quantum operation Ut- Note that we have implicitly subdivided the physical Hamiltonian into a term that 
cannot be varied, Hd, and a term that may be controlled externally. He- 
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FIG. 1: Bloch sphere representation of a single-qubit. 

To solve this control problem, we first transform Eq. into the following standard form: 

iU = {aa^ + vHc)U. 

This transformation is obtained by the following argument. Without loss of generality, we can assume that Hd 
o-iCx + a2(Ty + aaCz, where ax, cry, and cr^ are Pauli matrices. When + a| 7^ 0, let 

/ ai — a2« — fli + 021 \ 
- as /a + a^ 



(2) 
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2a 
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where a — y^jf-Fa| + af. It is straightforward to verify that fcfct — i.e., fc is a single-qubit operation. One can 
also verify that kHak^^ = auz- Then, letting Ui — kU, we have 

iUi = k{Ha + vH^)k-^Ui 
= (ao-^ + vkHck^^)Ui. 

where we have recognized that kHck^^ — H'^, which puts this into the form of Eq. The original control problem 
of steering Eq. from U{0) — I to Ut is thereby transformed into a problem of steering Eq. ||2Jl from Ui{Q) = fc to 
kUr- For control of a single qubit, we therefore only need to consider the standard form (0). 



A. Steering on Bloch sphere 



If no control signals are applied, i.e., v = 0, Eq. ([SJ simplifies to 

iU ^ aa,U. (3) 

The solution to Eq. Q is simply the rotation operator Rz{2at) — e~'°°'^*. We can thus easily implement any CTz 
rotation by turning off the control v and letting the system evolve for a certain time duration under Hd alone. We 
now define the notion of i?z-equivalence, namely, that Ui and U2 are i?z -equivalent if they satisfy 



(4) 



It is evident that this relation is reflexive, symmetric, and transitive, and therefore it is an equivalence relation on 
the Lie group SU(2) of all the single-qubit operations. To generate an arbitrary single-qubit operation Ut, it is then 
sufficient to generate an operation Ui that differs from Ut by a cr^ rotation from the right. Now our task has been 
reduced to a control problem on the space of all the i?^ -equivalence classes, that is, on a quotient space of SU(2). 

We now show that this quotient space is nothing but the Bloch sphere. Recall that the state of a single qubit can 
be represented by lip) = zi\0) + Z2\l), where Zi are complex numbers with restriction -|- |z2p = 1- Ignoring the 
global phase, this state can be written as 



cos-|0)+e^^sin-|l). 
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where the real numbers and define a point (x,?/, z) = {sin 6 cos (j>, sin 9 sin (j), cos 9) on a unit sphere, which is the 
well-known Block sphere shown in Fig. ^ To find and (f> explicitly, we can use the Hopf fibration, which is a map 
TT : 5*3 ^ 52 defined by = (x, y, z) e 5*2 16], with 

(z, x + yi)^ (|zip - |z2p, 2ziZ2). (5) 

We define a map $ fi'om a single-qubit operation to a point on the Bloch sphere by <&(C/) — U\l), where |1) is the 
single-qubit state corresponding to the south pole on the Bloch sphere. Then 

Furthermore, if Ui and U2 are i?^ -equivalent, we have $(C/i) = ^{U2)- Therefore, all operations that are j?^-equiyalent 
are mapped to the same point on the Bloch sphere. Conversely, it can also be proved by Lie group theory that 
all operations that are mapped to the same point on the Bloch sphere differ only by a Uz rotation from the right, that 
is, they are all i?z-equivalent. We thereby obtain that the quotient space of i?2-equivalence classes on SU(2) can be 
represented by the Bloch sphere. 

Now our control task becomes a steering problem on the Bloch sphere. Under the above mapping, the initial point 
of this steering problem is $(/) = and the final point is C/t|1), where Ut is the desired target operation. Our 
strategy is to first find the control v that can generate a trajectory on the Bloch sphere reaching the final point C/t|1) 
from the initial point Applying this control to the single-qubit system jSJ, we can obtain an operation Ut that 
is i?2-equivalent to the target operation Ut- The second step is to find the appropriate additional time duration tz 
such that Ut — UTe~^°''^'^' . The detailed solutions will depend on the form of the control term vHc- 



B. Oscillating control field 



We illustrate this approach of control on the quotient space with analysis of generation of single-qubit operations 
on two-level quantum systems subject to an oscillating electromagnetic control field A/2cos{Lot + 5). Some of the 
results shown below are known within the NMR literature, in particular that for a control field perpendicular to the 
static field. However, the solutions look different because the analysis here is made in the laboratory frame, instead of 
in a rotating frame as is customary in the NMR community. In NMR, under usual operating conditions the coupling 
terms simplify to an Ising interaction which is invariant under the rotating frame transformation. In the general 
case, when the coupling term is not of the Ising form, a time independent coupling Hamiltonian in the laboratory 
frame will be time dependent in the rotating frame. As noted above, implementation of a two-qubit local gate, i.e., 
ki = fell ® ki2, requires simultaneous implementation of each single-qubit operation fcn and ki2. We show below that 
such simultaneous single-qubit operations are readily steered in the laboratory frame. The same argument applies to 
local gates in larger multi-qubit arrays. 



1. Perpendicular control field 

We first consider the case in which the control Hamiltonian He is an oscillating control field perpendicular to a static 
drift field Hd. This is prevalent in many quantum systems, e.g.. Electron Spin Resonance (ESR). Here a single atom 
is placed in a uniform constant magnetic field directed along the z-axis, and the control is an oscillating magnetic 
field oriented along the x-axis [T^ . The dynamics of the quantum spin is determined by: 

tU - ya, + ^ cos{iut + j U. (6) 

A solution to Eq. (01 for near resonance, u; ~ wg, may be found by transforming first to a rotating frame defined 
by the drift (Larmor) frequency, i.e., Ui = g,'^^ot/2a^jj ^ -y^g then obtain 

A 

iUi = — cos{LLjt + 5)(^ax cosujot — (jy sinujQt^Ui 

A ( gi[(w+wo)t+<5] ^ g-i[(w-wo)t+<5] 

~ ~^ I g-i[(w+(^o)t+'5] _|_ gi[(i^-i^o)*+'5] 

Since for near resonance between drift and control fields the terms ^^A{^+^oY+S\ oscillate much faster than the terms 
g±i[(w-wo)t+5]^ they make on average little contribution to U\ and can be neglected (rotating wave approximation). 
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FIG. 2: The effect of the single-qubit operation U(tf,tz), Eq. IIOL on the quantum state |1): (A) u 7^ uoo', (B) to = ljo- 



resulting in the approximate equation: 



iUi = — [ cos[(aj — i^o)^ + S\ax + sin[(ti; — i^o)^ + S\ay ) Ui. (7) 



Now letting U2 = e»[("-'^o)*+'51/2'^- [/i, we have 



«C/2- -c7, + ^— a, C/2, (8) 



which has an expUcit solution 

U2{t) = g-i(A/A<J^ + {uo-u)/2a^)t ^ 

Combining these two transformations, we obtain an approximate solution (i.e., within the rotating wave approxima- 
tion) to Eq. lO as 

U{tf) — g-i('^*/+'5)/2f^g-»('4/4(T^ + (cJo-i-^)/2(T^)t/gii5/2(T3 ^ ^g-j 

This is equivalent to the standard treatment of ESR and NMR • Since we can also generate any az rotation by 
turning off the control field for a time duration t^, a general form of the quantum operation generated by this 
system can then be written as 

U{tf,tz) = C/(t/)e~*"''/^'^^*^ = e~*("*^+^)/2'^^e~*^^/'''^=""'"^"''~")/^'^'^*^e*^^~"°*^^/^'^' . (10) 

In the language of control on the quotient space developed above, the effect of the operation U{tf,tz) on the 
quantum state |1) can be described directly from Eq. H10|l as a rotation about the axis {A/2,0,ujq — uj) followed by 
another rotation about the z-axis, as shown in Fig.|21for both near resonant (A) and resonant (B) situations. From 
this it is evident that when lu ^ loq, all points (x, y, z) on the Bloch sphere satisfying 

shown as the shaded area in Fig. [2 can never be reached for any choice of tuning parameters A, 5 and time durations 
tf, tz- Consequently, the reachable set in this near-resonant case is not all single-qubit operations. To allow imple- 
mentation of any arbitrary single-qubit operation, we must impose the strict condition uj — coq, that is, the oscillating 
field must be tuned to the resonant frequency. 

We now show that when operating at exact resonance, the time duration to implement an arbitrary single-qubit 
rotation in the laboratory frame is quantized. At resonance, Eq. H10(l becomes 

U{tf,tz) = g-i(^otf+S)/2a^ ^-iA/4a^tf ^i{S-uot^)/2a^ ^ j--|^2) 
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On the other hand, we know from Euler's ZXZ decomposition that an arbitrary single-qubit operation Ut can be 
decomposed as 

Note that in Eq. H13() the coefRcients in the first two terms are chosen such that Ut is mapped to the point {9, 4>) on 
the Bloch sphere. Comparing Eqs. (|f 2|l and we find that for oj — wq, an arbitrary single-qubit operation Ut in 
Eq. 1)13(1 can be implemented if the following conditions are satisfied: 

(/) - 7r/2 + 2mi7r = wot/ + (5, (14) 

7r-6' + 2m27r = -tf, (15) 

-7 + 2m37r = (J-cJo^^, (16) 

where mi, m2, e Z, and mi + m2 + ms is even. Any values of the control parameters A and 5, and time durations 
tz and tf that satisfy Eqs. (|14|l - (|16() can generate the desired target operation Ut within the total time tf +tz. From 
Eqs. (dl) and (QH)), we also find that 

+ 7 - 7r/2 + 2m7r 
tf+tz = , (17) 

UJQ 

where m is an integer. Hence, from Eq. (|17|l . it is evident that if the oscillating field is tuned at resonant frequency, 
the time duration to implement a single-qubit operation is quantized with a time duration that is independent of 
the amplitude of the oscillating control field. This is a particularly important issue in the generation of two-qubit 
local unitary operations, which are given by ki — kn (g) fci2, where fcn and fci2 are both single-qubit operations. To 
implement such a two-qubit local unitary fci, we need to implement both fcn and fci2 within the same time. We 
can satisfy this constraint by making the oscillating control field that addresses one qubit to be resonant and that 
which addresses the other qubit to be slightly off-resonance. Using this offset as an additional variable, we can then 
apply a numerical optimization procedure to achieve fcn and fci2 simultaneously, leading to the desired local two-qubit 
unitary. 

The above analysis of the time duration is made in the laboratory frame and thus differs from the standard 
NMR/ESR situation in which pulses are made in the rotating frame 18J (Eq. 1)17(1 applies only in the laboratory 
frame). When working in the rotating frame Ui = ^li^ot / 2<j ^ -jj ^ usual in NMR experiments at resonance, the 

following time constraints are instead required to implement Ut'- 



In this case, we have 

tf+tz 



0-7r/2 + 2mi7r = 5, (18) 

Ti-e + 2m2TT = -tf, (19) 

-7 + 2TO37r = 8-u}Qtz. (20) 

2(7r -9) + 4m27r + 7 - 7r/2 + 2m7r 
A Wo 



Hence, when working in the rotating frame, it is possible to decrease the time to reach a desired operation merely by 
increasing the amplitude of the oscillating field. 



2. Non-perpendicular control field 



We now consider the situation in which the control Hamiltonian is an oscillating electromagnetic field that is not 
perpendicular to the static field: 



iUi 



LOq 



A 



cos{iLjt + S){cos(^az — sin^CTa;) 



Ui. 



(21) 



Here the control Hamiltonian He = cos (az — sin Cctj;, so the static and oscillating fields are tilted by an angle. This 
can arise in NMR or in coupled semiconductor quantum dots jl9(|. Let U2 = e^'^'^o^+^y^'^'^Ui. Imposing the rotating 
wave approximation again, we have 



A 

iU2 = — [2 cos ( coa{ujt + 6)az — sin C(cos(ti; — uJo)t + sin(a; — iUo)t CTy )] U2 



(22) 
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When the oscillating field u is tuned at the resonant frequency wq, Eq. 1)22(1 becomes 

A 



iU2 — — [2 cos Ccos{uJot + 6)az — sm(ax]U2. (23) 



Letting t/3 = e-"/4'^''t/2, we transform Eq. (|23|) into the standard form: 



A 

iUy, = — \2 cos C cos(a;o^ + 5)ax + sin Qaz\ U3. (24) 



Finally, letting U4 = e'^sinC/4<T,t[/g^ obtain: 



v4 

if74 — — cos C cos(tJo* + 



o'x cos ( — sin (t) — Gy sin ( — sin C,i) 



U4 



A f gi[(a;o-A/2sinC)t+<5] _|_ gi[(-aJo-A/2smC)t-<5] 

JCOSC Li[(a;o+A/2sinC)t+<5] _,_g' 



V'^O^C 1 „i[(wo+A/2sinC)t+<5] I „i[(-wo+A/2sinC)t-<5] ) ^4 



(25) 



Since the control fields are usually high frequency signals, we have ujq 3> A/2 sin^. Eq. (|25|l can then be simplified to 

A 

iil^ = — cos (cos{ujot + S) (TxU4, (26) 

which has a solution 

U4{t) = g^'-4 cos C[sin(LJot+<5)-sin i5]/(2wo)o'x (27) 

Combining the above series of transformations, we obtain an approximate solution to Eq. H21|l as: 

Uitf) — g^*('^ot/+'5)/2(T^giAsinC/4(Ta:t/g-i[AcosC[sin(uJot/+'5)-sin(5]/(2(^o)-'5/2]crj ^28) 

As in the previous example, we also have the freedom to implement any arbitrary cr^ rotation by turning off the 
oscillating control field for a time duration tz. Hence the final quantum operation achieved can be written as: 

U{tf tz) = U{tf)e~^^''^'^'^'*' = g-*('^ot/+'5)/2o-2giAsinC/4o-^t/g-i[AcosC[sin(cJot/+<5)-smi5]/(2i^o)-'5/2+i^ot2/2]tTj ^29) 

The time duration required to implement an arbitrary single-qubit operation under these conditions may be similarly 
derived by comparing with the general form of Ut given in Eq. \i'6\ . This results in the conditions: 

(/) - 7r/2 + 2mi7r = ujctf + S, (30) 

■7T — 9 + 2m2TT = ——sm(tf, (31) 

— 7 + 2m37r = -[smiujotf + S) — sinS] + S — uatz, (32) 

where mi, m2, G Z, and mi + m2 + ma is even. From Eqs. H3l)|l and (|32|l . we then find that 



tf+tz 



TT ACOSC, , . c.\ 

+ 7 1 (cos0 + smd) + 2m7r 

2 ujQ 



LOo, (33) 



where m is an integer. Thus when the oscillating control field is tuned at the resonant frequency, the total time to 
implement a single-qubit operation is also quantized in this case of a non-perpendicular control field. 



C. Constant control field 



We now consider the case when the control v can be varied within a given range [wojWi]. A well known control 
strategy is Bang-Bang control, which pertains when the control field switches back and forth between two extremal 
values vq and vi • In classical control theory, it can be shown that Bang-Bang control is indeed the time optimal 
control strategy for a double integrator system. We will show that for a quantum system of the form in Eq. 
Bang-Bang control can generate an arbitrary single-qubit operation with the minimum number of switchings. We 
note that the term Bang-Bang control has recently been adopted with a somewhat different meaning in the study of 
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dynamical coupling of open quantum systems [20I I21I , where it is referred to performing instantaneously or 

as fast as physically possible, a set of unitaries, implying that the corresponding set of Hamiltonians can be turned 
on for negligible amounts of time with (ideally) arbitrarily large strength. A modified form of this approach for finite 
control field amplitudes and durations has been presented in Ref. p^ . 
We rewrite the quantum system (0) as follows: 

iU = {aa, + vHc)U, (34) 

where v € [vo,Vi]. Without loss of generality, we can assume that vq = 0. By taking extremal values v = and 
V = vi, we obtain two linearly independent Hamiltonians 

Hi = aaz, (35) 
H2 = aaz + viHc- (36) 

Our task is to achieve an arbitrary single-qubit target operation Ut by switching between Hi and H2. The final 
operation generated by switching between these two Hamiltonians can be described as 

In the mathematical literature, this problem is also known as the uniform finite generation problem on Lie groups, and 
can be traced back to the early 1970's H^H^- Recently, D'Alessandro studied the optimal evaluation of generalized 
Euler angles on the rotation group with this approach '27] . 
To solve this problem, we first rewrite H2 as: 

H2 = aoz + viHc = hia.j, + 620'^ + ^z^z- 

Let k — e*'^^^, where 7 satisfies 

-hi sin 27 + &2 cos 27 = 0. (38) 
From the Campbell-Baker-Hausdorff (CBH) formula we have: 

Hi = kaOyk'^ — acTz, 

H2 = kH2k^ = e'^'-'ibia^ + &2a.y + b3az)e-'^'^ 



= (61 cos 27 + 62 sin 27)CTa; + b^az 
~ 6(sina(T^ + cos aax), 



(39) 



where 



and 



6= -^(61 cos 27 + &2 sin 27)2 + 6§ (40) 



{Hi,H2) tr{HiH2) tT{HiH2) 
sma = — : — : — r = = 5 = r = 1 r- 41 

(Hi,Hi)^H2,H2)i tr(i/2)^tr(i/|)5 tr(i/2)^ tr(H|)5 

Without loss of generality, we can take a E [— ^, When the target operation Ut can be generated in Bang-Bang 
strategy as in Eq. H37() . we have 

kUrk^ = ^e~^ffl^„g-^ff2t„-l ...g-»i/it3g-«ff2t2g-iffiti;.t 

_ gi//it„giH2t„^l . . . giHit'i ^iH2t2 ^iHiti 

In order to generate the target quantum operation Ut from the original quantum system H34|l . we therefore only need 
to generate kUTk^ from the two Hamiltonians Hi and H2 in Eq. (|39|l . For the special case a = -I, we have H2 = ba^ 
and Ut simplifies to the familiar Euler ZXZ decomposition of Eq. H13|) : 

The Bloch sphere steering technique can now be used to find the time sequences for switching between Hi and H2 ■ 
We know that the effect of e~*^i* on the quantum state |1) is a rotation about the z-axis, whereas the effect of e~*^^* 
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FIG. 3: Effect of the rotations e~'^^^ (A) and e~'^^* (B) on the state of a single qubit. (A) a rotation about the z-axis; (B) 
a rotation about the axis n — (cos ^, 0, sin (C) successive switching between the two rotations (A) and (B) applied with 
finite time durations allows any point on the Bloch sphere to be reached. 



is a rotation about the axis (sin a, 0, cos a), as shown in Fig. |3[A) and (B). Application of the time sequence given in 
Eq. (|37|l to the quantum state |1) results in an alternation of rotations about these two axes, illustrated in Fig.|2tC). 
Therefore, to generate Ut, we need to reach the point Ut\1) from the south pole |1) on the Bloch sphere by switching 
back and forth between two rotations, ft is evident that we can now achieve any arbitrary point on the Bloch sphere, 
by using these two rotations in alternation, with finite switching times. Solution for the number of switches and time 
durations is illustrated explicitly by the following example. 

Example 1 (Bang-Bang solution for Hadamard gate) 

Suppose that in the single-qubit system (0), we have V e [0,a\/3]. This gives the following two Hamiltonians defined 
by the extremal values of v: 

Hi = aaz, 

, TT , (42) 

Hi = 2a(sin -a^ + cos -da;). 

6 6 



We aim to generate the Hadamard gate 



by switching between these two finite Hamiltonians at well defined times. We first find the point corresponding to 
the Hadamard gate Uh on the Bloch sphere. Computing 

and using Eq. (0, we find that Uh corresponds to the point (0,0) — (§,7r) on the Bloch sphere. The Hamiltonians 
Hi and i?2 of Eq. (|42|l generate two rotations on the Bloch sphere, one about the z-axis and the other about the axis 
n = (cos ^, 0, sin ^) (see Fig.0J. We can implement the Hadamard gate Uh by switching between these two rotations 
as follows: 

Uh = -ie-^"^^^e-'"-'^^e-'"^^\ 

where ti = 9.90/a and t2 ~ H — cos^^ From Fig. 01 it is also clear that for this Hamiltonian, we can generate 
any arbitrary single-qubit operation with at most three switchings. 

The minimum number of switchings between the two Hamiltonians Hi and H2 needed to generate a target quantum 
operation Ut is determined by the rotation axis (sin a, 0, cos a) and the z coordinate of the point C/t|1) on the Bloch 
sphere. When a = 0, we can get to any point on the sphere starting from the south pole (|1)) by at most two 
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FIG. 4: Two rotations generated by the Hamiltonians given in Eq. I42II . 

switchings. This corresponds to the Euler ZXZ decomposition. For a general value of a € [0,7r/2), it is not hard 
to derive that the coordinates for A„ in Fig. OIC) is A„ = [cosq;„, — sina„], where q;„ — 11/2 + 2{n — l)(7r/2 — a). 
To ensure any point on the sphere can be reached by the combination of these two rotations, we only require the 
condition that a„ > 37r/2, which in turn yields that we can implement any arbitrary single-qubit operation with at 
most [7r/(^ — a)~\ switchings, where the ceiling function \x~\ is defined as a function that rounds x to the nearest 
integer towards infinity. Similarly, when a € (7r/2, tt], the maximal switchings needed is [tt/ (a — ^)] . To achieve the 
minimum number of switchings, we therefore need to use the extremal values of the control to minimize the angle a. 

We note that the generation of arbitrary single-qubit operations by switching between two constant control fields 
has many other applications in quantum computation and quantum information processing. For example, in quantum 
simulations, if we are given two single-qubit Hamiltonians Hi and H2, the usual way to generate U = e'^'^* is to use 
the Trotter formula plIsollsH 

We can use the above Bang-Bang approach to simulate U — e*^^* in a constructive and exact manner, without 
assumption of impulsive and arbitrarily strong controls 2Qj| or of infinitesimal operations [33 | . Another example is in 
encoded universal quantum computation with the exchange interaction [s^l ■ Here the encoded single-qubit quantum 
gates are generated by two Hamiltonians identical to those in Eq. H42|l . The switching sequence to implement any 
such encoded single-qubit operation can readily be derived using the strategy outlined above, which yields the result 
that only four switches are necessary, corresponding to the solution found by numerical arguments in js^. Finally, 
it has been shown recently that electron spin in a semiconductor heterostructure may be controlled by a single-qubit 
Hamiltonian generated by g-tensor engineering that is realized by applying a DC bias voltage 34J. By switching 
between two Hamiltonians obtained from two different g- factors corresponding to different bias voltages, we can thus 
implement any single-qubit quantum operation. 

III. GENERATION OF TWO-QUBIT OPERATIONS: STEERING IN THE WEYL CHAMBER 

To generate an arbitrary quantum operation on arbitrarily many qubits, the target quantum operation can be 
decomposed into a combination of single- and two-qubit operations 1]. In this section we now consider the generation 
of two-qubit operations. 

A general Hamiltonian for a two-qubit quantum physical system can be written as 

H^gi-~a(»I + I(^g2-~^ + S, (43) 

where 'a is the vector {ax, <Jy, <Jz) of Pauli matrices, and S = ^ Jai3<j\a^p with (j]^(J^p = aa crp, a, P G {x, y, z}, and 
Ja(3 the corresponding coupling strength. The coefficients gi = [gix, 9iy, 9iz], 92 = [92x, 92y, 92z] may be regarded as 
single-qubit "g- factors" . alt houg h they may have very different physical origins, depending on the physical realization 
of the qubits. From [lOl llsL l29l | . we know that for any arbitrary set of coupling strengths {Ja/d}, there exists a local 
operation k such that 



kSk^ = Jxcrlal + Jycrlal + Jzala^. 



(44) 
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FIG. 5: Tetrahedral (or Weyl chamber) representation of local equivalence classes of nonlocal two-qubit operations from 
Ref. [T^. Points within the tetrahedron are labeled [ci, 02,03], where tt — C2 > ci > C2 > C3 > 0. Points O([0,0,0]) and 
Ai{[tv,0,0]) correspond to local operations, L([|,0,0]) to the CNOT gate, and Ag^^, f , §]) to the SWAP gate. 

Without loss of generality, we can assume for the two-qubit steering problem that the two-body interaction term S is 
always given as in Eq. (|44|l . We shall study the generation of an arbitrary quantum operation from relevant physical 
Hamiltonians by tuning of the parameters gi, g2, and Jx, Jy, Jz, all of which correspond to control parameters in this 
context. 

In Ref. [l3| we have shown that this problem can be reduced to a steering problem in a tetrahedron. Central to 
this analysis is the notion of local equivalence. Two quantum operations U, Ui £ SU(4) are called locally equivalent if 
they differ only by local operations: U = kiUik2, where ki, k2 G SU(2) ® SU(2) are both local operations. Clearly this 
defines an equivalence relation on the Lie group SU(4) of all the two-qubit operations. It was shown in ,15] that the 
local equivalence classes of two-qubit operations are in one-to-one correspondence with the points in the tetrahedron 
OA1A2A3 as shown in Fig. |S1 except on its base (where there is a one-to-two correspondence). This tetrahedron 
is a Weyl chamber that contains all possible local equivalence classes. To locate the point in the tetrahedron that 
corresponds to a given two-qubit operation, we make use of a set of local invariants. Ref. provides a simple 
procedure to compute three real quantities gi, (72 and gs for a given two-qubit operation and showed that two 
quantum operations are locally equivalent if and only if they have identical values of these three invariants. These 
local invariants can be expressed in terms of the tetrahedron coordinates [ci , C2 , C3] as follows jl5| | : 

gi — 4: cos Cl cos C2 cos C3, 

(72 = 4sinci sinc2 sinc3, . (45) 

33 = cos 2ci -I- cos 2c2 + cos 2C3 

By solving Eqs. I|45|) . we can obtain the point [ci, C2, C3] in the tetrahedron that represents a given two-qubit operation. 

As time evolves, the quantum system controlled by H generates a continuous flow in the space of all the quantum 
operations. At any time instant t, the quantum operation U{t) on this flow can be mapped to a point in the 
tetrahedron. This defines a continuous trajectory within the Weyl chamber. We then employ a two step procedure 
to generate a target two-qubit operation from a given Hamiltonian. First, we steer the Weyl chamber trajectory 
generated by the given Hamiltonian H to reach the point representing the local equivalence class of the desired target 
operation. Second, we perform the necessary local operations that transform the quantum operation at this Weyl 
chamber point to the true target operation. The single-qubit operations are generated as described in Section ^] 
Therefore we only deal with the first step here, namely steering of the system Hamiltonian H to generate a trajectory 
that hits the target in the tetrahedron OA1A2A3. 

For a general Hamiltonian containing both local and nonlocal terms, it is usually difficult to implement an arbitrary 
two-qubit operation directly. This is reflected in the complexity of a general trajectory in the Weyl chamber. One 
approach is to introduce an intermediate step to construct a universal gate set of elementary gates that can be used 
to construct any arbitrary two-qubit operation 0. An alternative perspective is to seek to minimize the time for 
nonlocal operations |10l l35| . We follow the first approach of steering to achieve a universal gate here. The solutions 
will then not necessarily be time optimal. Optimization for time can be performed subsequently in comparison of 
different gates |36j . 

One well-known example of such a universal gate set is the CNOT gate together with single-qubit operations 0]. 
In [33, we showed that both CNOT and Double-CNOT (or its locally equivalent variant, iSWAP 38, 39]) require at 
most three applications to realize any two-qubit operation. More recently, we discovered a new quantum operation 
B that is universal by at most two applications together with at most six single-qubit operations [361 ]. Whenever the 
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direct generation of any quantum operation is difficult, we can take advantage of these results for efficient two-qubit 
quantum circuits and first seek to implement the quantum gate B or CNOT, then use the corresponding analytic 
circuits of |2Ein3 in order to construct an arbitrary two-qubit operation. 

To illustrate the basic idea of the steering approach, consider first a simple case when gi — g2 — 0, i.e., the 
Hamiltonian contains only the nonlocal term S. The Weyl chamber trajectory generated by this Hamiltonian is the 
straight line U{t) = [Jx, Jy, Jz]t p^ . The evolution direction of this line can be changed by conjugating with an 
appropriate local operation from the Weyl group. Since the tetrahedron is a three dimensional geometric object, it is 
evident that by changing the direction at most twice, we can steer the Hamiltonian to anywhere in the tetrahedron. In 
other words, by turning on the Hamiltonian at most three times, together with application of at most four local gates, 
we can implement any arbitrary two-qubit operation [Tsl l . Note that when the local terms commute with the nonlocal 
term, we can drop the local terms and the trajectory is unchanged. For example, with H — gixCrl + g2x'^x + Jx^'x^^'^, 
corresponding to the Hamiltonian derived in the superconducting qubit proposal of ^10], the Weyl chamber trajectory 
is exactly the same as that generated by JxCr^a^, namely along the straight line OAi. In this case, the straight line 
trajectory leads directly from the origin (local unitary) to the CNOT gate ,0,0]. This result was derived by more 
complex means in Ref. . 



A. Isotropic qubit coupling 

Isotropic coupling of qubits, characterized by coupling amplitudes Jx = Jy = Jz ^ J is a, very common form of two- 
body interaction that appears in many physical qubit proposals involving electronic spins, e.g., spin-coupled coupled 
quantum dots t41=)42| and donor spins in semiconductors 43, 44, 45, 46, 47] . We present here two different strategies 
to steer the Weyl chamber trajectories generated by this Hamiltonian, both of which are based on the ability to tune 
the relative values of the single-qubit parameters gi and q^. The technological ability to make this tuning has been 
demonstrated for electron spins in semiconductors in Ref. j34l ] . For greater simplicity, in the remainder of this section 
we use a slightly modified notation for the two-qubit interaction, namely the form JS where S — cr^.a'^ + Uyay + a^a^ 
and J is the isotropic coupling amplitude. 

The first strategy is based on tuning gi and g2 to be equal, i.e., setting gi — g2- Then the individual components 
of the total single-qubit spins commute with the interaction term S, 

\a\ -f ul, S\ = 0, {a\ + al, S\ = 0, \a\ + al,S\ = 0. 

This implies that single-qubit and two-qubit terms commute, 

[51 • "(T (g) / + / gi • "(T , 5] 
= \9\x{(y\ + crl) + giy{al + al) + 5i^(ct] + (^l),S] = 0, 

allowing the time evolution to be factorized exactly: 

^iHt _ ^igi-~at ^ ^ig2-~at , ^iJ St ^ 

Consequently, the flow generated by the Hamiltonian H in the tetrahedron is the same as that generated by S, i.e., 
we have achieved a reduction to the purely nonlocal interaction. The resulting Weyl chamber trajectory lies along 
the line OA3A1, and we can achieve the CNOT gate by changing the evolution direction just once (cf. Example 4 
in Ref. ^3)- emphasize that this reduction results from the tuning to equivalence of the two single qubit terms. 
Furthermore, it is valid independent of the nature of the single-qubit terms, i.e., whether they involve ax, cry and/or 
az- Thus, although the result resembles the known result that the isotropic exchange interaction generates a vSWAP 
gate from which a CNOT gate can be obtained by conjugating with single-qubit gates ^i§\, it is valid here for a more 
general Hamiltonian that contains both nonlocal and local terms. The extended range of validity is a consequence of 
the decomposition of the two-qubit operation into nonlocal equivalence classes and local operations. 

The second strategy is based on tuning the single-qubit parameters gi and g2 such that they become linearly de- 
pendent, i.e., gi = Xg2, where A describes the parameter ratio. This is clearly a less stringent tuning requirement than 
the first strategy. After some mathematical analysis, the trajectory generated by the Hamiltonian in the tetrahedron 
is found to be: 

ci = 2 Jt, 

2 J 

C2 = — \ sin^ ( — sinwi)!. 
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FIG. 6: Weyl chamber trajectory obtained upon steering the isotropic Hamiltonian to reach CNOT by tuning the local terms 
gi and g2 to be linearly dependent. The Weyl chamber trajectory is located in the OA3A1 plane of the tetrahedron, and goes 
from the origin to the point [7r/2, 0, 0] representing the local equivalence class of CNOT, with a total of four oscillations. 



where w = v/(A - l)2||g2|p + 4J2. Since C2 = C3, the entire trajectory lies in the OA3A1 plane of the tetrahedron 
OA1A2A3, as shown in Fig.lHI CNOT is located on the intersection of this plane with the base of the tetrahedron, at 
[tt/2, 0, 0], and hence constitutes a natural target for these trajectories. It is readily verified that in order to steer the 
trajectory to CNOT we only need to satisfy the conditions 

2Ji = |, (46) 
ujt ~ rmr, (47) 

where m is a positive integer. From these two equations, we find that the time required to reach CNOT with this 
trajectory is t — jj, and the parameter ratio A has to satisfy the condition 

(A-l)'||.g2||' = (16m2-4)J^ (48) 

From Eq. H46I) . we find that the time to achieve CNOT depends only on the coupling strength J, while Eq. I|47|) 
shows that the integer m determines the number of oscillations in the trajectory. It is interesting to note that this 
time is equal to the time needed to achieve CNOT from a purely isotropic Hamiltonian JS, or from the first strategy 
discussed above, both of which have a straight line Weyl chamber trajectory [l5|. 

As an example, consider a system with J = 0.1 and choose m = 4. Then Eq. H48fl becomes 

(A-1)2||52|P = 2.52, 

which needs to be solved for A, given a specific set of g parameters. Choosing 172 = [4,4,4], we obtain A = 0.7709, 
uj — 1.6, and find a time t = 2.57r to reach CNOT. The resulting Weyl chamber trajectory is shown in Fig. 



B. Anisotropic qubit coupling 

Hamiltonians of the general form of Eq. I|43|l with anisotropic coupling coefficients Jap are encountered in many 
proposed solid state implementations of quantum computation. These include quantum dots |49j| , electrons coupled 
by long range quantum Hall effects , electrons on helium 'si'l , and atoms in cavities [s^l • 

To illustrate the Weyl chamber steering approach for anisotropic qubit couplings, we consider here a Hamiltonian 
containing an Ising interaction, i.e., a single diagonal component of Jap only: 

Hyy = gixcrl + gizc^l + g2xcrl + g2z(jl + J^yl^yl- 

This Ising-type coupling between qubits is commonly seen in proposals for superconducting Joscphson junction 
qubits 0, I53I and also arises as one limit of dipole-dipole or J-coupling systems 0, 0|. It is straightforward 
to show that the Hamiltonian Hyy is locally equivalent to the corresponding Hamiltonians H^x and H^z '■ 

Hxx = giyC^l + giz<jl + g2ycrl + g2zcrl + Jo-lal 

Hzz = gix<ji + giy^l + g2xol + g2y(7l + Jo-\(jI. 
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FIG. 7: Minimum time Weyl chamber trajectory that reaches the B gate from the Ising-type YY interaction in one switching. 
See Fig. for definitions of the coordinates and special points. 



Consequently the Weyl chamber trajectories generated by H^x and H^z are the same as those generated by Hyy. 
Using Makhlin's procedure |3 to evaluate the local invariants of U{t) — e'^""*, we find 

_ {fl + J^)J^x^ + (/I + .P)JW + /f/l - J' 

92 = 0, (49) 
53 = (s/f /I - J'if? + fi) + J\^^W + 3) + ^JW 

■ (/I - J2) + U'x\fl - J2)) / [fl + J2)(/| + J% 

where x = cos ^/ff+J^ t, y ^ cos V/i + t, fi = a/sL + ffL + VfL and /a = + " V^L 

Without loss of generality, we can assume that J = 1 in the subsequent analysis. 

It turns out that both CNOT and B gates can be generated by turning on this anisotropic Hamiltonian only once. 
The B gate is more efficient than the CNOT gate in representing arbitrary two-qubit operations 36] and may therefore 
be more useful for quantum simulations, while the CNOT gate is required for many common error correction protocols. 
We first consider realization of the B gate. The local invariants are 5i = 32 = ff3 = [s^. To find the Weyl chamber 
trajectory that reaches the point B — [7r/2, 7r/4, 0] in time t we need therefore only to set gi — = in Eq. 
This yields the following equations for t that must be satisfied simultaneously: 

cos 2^/2 + 1 < = _/2±^(/2 + l), 

cos2^/| + l t = -fi T ^{fl + !)• 

There are infinitely many solutions to these two equations, depending on the combination of single-qubit parameters 
gia- Numerical analysis reveals that the time optimal solution for reaching the B gate is achieved by setting /i = 
1.6753, /2 = 0, with terminal time t = The corresponding Weyl chamber trajectory is shown in Fig. [7| Note that 
this time optimal trajectory remains at all times in the basal plane of the tetrahedron OA1A2A3. A similar analysis 
was made in |15| for a particular instance having gi^ = g2x — —0^/2, giz = g2z = 0, and J = a^. 

We next consider realization of the CNOT gate. The local invariants for the CNOT gate are gi — g2 ~ and 
gs — I 01 . Solving Eq. (|49|l we obtain the following two equations: 

COS2y^/2 + l i = 

cos2y^/| + l i = -fi. 

By numerical analysis, we find that the minimum time to achieve CNOT is t — 3.2551, with /i — 0.9516 and 
/2 = 0.9492. Comparing this with the time optimal solution for reaching the B gate described above, it is clear that 
for such Ising-type coupfing, the B gate can be generated in shorter time than the CNOT gate. 
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FIG. 8; Weyl chamber trajectory designed to reach CNOT and generated by the weak two-qubit interaction given in Eq. 1)52^ . 
The trajectory is at all times in or very close to the OA^Ai plane of the tetrahedron OA1A2A3 of Fig. |K| 



C. Weak qubit coupling 

In many coupled qubit systems, the two-qubit interaction term is much weaker than the local term, that is, J is 
small in comparison to gi and 32- For example, this is the case in circuits of flux qubits coupled by a mutual inductance 
that is controlled by the circulating current in a DC Superconducting QUantum Interference Device (SQUID) |0. It 
is also the case in most NMR systems of exchange-coupled nuclear spins . Estimates of parameters for exchange- 
coupled electron spins in quantum dots also put these systems in this weak coupling regime at magnetic field strengths 
relevant to experiments |42L Issl 15(1 Is^ . In such situations we have found that the Weyl chamber trajectory can 
be approximated to good accuracy by some simple curves within the tetrahedron. This provides a useful strategy 
in many physical systems, where instead of attempting to control the exact time evolution, we seek rather to steer 
only the approximating curves in the Weyl chamber and may wish to incorporate additional factors into numerical 
optimization protocols. 

We first consider a simple case, namely, tuning gi and 52 such that their norms are equal, i.e., \\gi\\ = \ \g2\\- In this 
case, when in addition J <C \ \gi\\, the trajectory generated in the tetrahedron can be well approximated by a straight 
line. Specifically, when the Hamiltonian is given in the Ising form 

H^gi-~a(E}I + I(E}g2-~a + Jz^lcrl 

with 1 1 (71 1 1 = II52II and Jz small, the resulting Weyl chamber trajectory can be approximated by the following straight 
lines: 

Jzt 



[C1,C2,C3] = 



\9i\? 



[2gizg2z,p,p], if 251^522 > p; 



m\\ 



where p — y gf^ + 9iyy 92x + 92y Recall that in the case when the Hamiltonian is purely nonlocal, i.e., it contains 
no singie-qubit terms, the Weyl chamber trajectory is a straight line and we can generate any arbitrary two-qubit 
operation by applying local unitaries to change the evolution direction of the trajectory. We can use the same technique 
here to steer the Weyl chamber trajectory to arrive at any point in the tetrahedron with at most two switchings. 

For the general case when \\g\\\ ^ II52II and S = Jxcr^crl + JyO'yfJy -\- Jzcr^cr^, we find that the trajectory can be 
well approximated by the following sinusoidal curve in the tetrahedron: 

_ Jxgixg2x + Jygiyg2y + Jzgizg2z ^ 

II51II -115211 ' (50) 

C2 C3 = p{gi,g2)\sm{\\gi\\ - \\g2\\)t\, 

where p is a nonzero function of gi and 92- From Eq. (|50|l . the approximate curve stays in the plane OA3A1 and moves 
away from the origin in the direction of CNOT. Therefore, for this Hamiltonian, the CNOT gate is a natural choice 
of target operation. To reach CNOT, we require that Eq. 1)50(1 be satisfied with values ci = ^ and C2 = C3 = 0. It is 
clear that the form of the function p is irrelevant. Consequently, we need only to satisfy the following two equations: 

Jxgixg2x + Jy9ly92y + J zgizg2z r,^ _ TT 

II51II-II52II "2' (51) 

(Il5i||-Il52||)i = m^, 
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where m is an integer. To illustrate this approach, consider the following example of weak coupling: 

H = gi-~a ® I + I ® g2-~o' + Q.2<tI(jI. (52) 

Here the local terms gi and 52 are one order larger than the interaction strength Jz — 0.2. From Eq. ()51|1 . we find 
that the local terms gi and .92 must satisfy 

II51II -113211 -(115111-115211) 4m- ^ ^ 

Setting, e.g., m = 3, we derive one solution for Eq. (|5!^ as 51 = [2.5,0,10.0182] and 52 — [2,0,7.8177]. After a time 
duration t = 4.1778, this Hamiltonian can achieve a gate located at [0.5000, 0.0002, 0.0002]7r, which is very close to 
the CNOT gate. The Weyl chamber trajectory is shown in Fig.|Sl where the solid line is the real trajectory generated 
by the Hamiltonian, and the dotted line is the approximate sinusoidal curve. A realistic application of this example 
was recently made in Ref. |1^ to inductively coupled flux qubits for which ||gi|| ||g2|| and Jx = Jy = 0. 



IV. CONCLUSIONS 



We have developed a systematic approach to generate arbitrary two-qubit quantum operations by steering the 
unitary evolution corresponding to local and nonlocal interactions. The generation of single-qubit quantum gates was 
shown to map onto a steering problem on the Bloch sphere. Local two-qubit gates corresponding to coupled single- 
qubit gates can then be generated by combined steering on two coupled Bloch spheres. The nonlocal components of 
two-qubit gates were shown to be analyzable as a steering problem in a Weyl chamber of local equivalence classes. 
We applied this steering approach here to the generation of two-qubit quantum gates in several physical examples 
relevant to solid state implementations of quantum information processing. The methods described here are expected 
to provide useful guidance to experimental design of circuits and pulse sequences for realization of quantum logic 
gates. 
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